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ABSTRACT 



Context. Previous studies of type IIP supernovae have inferred that progenitor masses recovered from hydrodynamic models are higher than 
Q 15 M Q . 

S— j Aims. To verify the progenitor mass of this supernova category, we attempt a parameter determination of the well-observed luminous type IIP 
r/3 | supernova 2004et. 

C3 , Methods. We model the bolometric light curve and the photospheric velocities of SN 2004et by means of hydrodynamic simulations in an 
extended parameter space. 

I , Results. From hydrodynamic simulations and observational data, we infer a presupernova radius of 1500± 140^q> an ejecta mass of 24.5 + IMq, 
^» ■ an explosion energy of (2.3 ± 0.3) x 10 51 erg, and a radioactive 56 Ni mass of 0.068 ± 0.009Mq. The estimated progenitor mass on the main 
£f) ' sequence is in the range of 25 - 29Mq. In addition, we find clear signatures of the explosion asymmetry in the nebular spectra of SN 2004et. 
f*"*) , Conclusions. The measured progenitor mass of SN 2004et is significantly higher than the progenitor mass suggested by the pre-explosion 
■ images. We speculate that the mass inferred from hydrodynamic modeling is overestimated and crucial missing factors are multi-dimensional 
CN ' effects. 
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1. Introduction light curve and the evolution in the photospheric velocities re- 
quires both high-quality photometric and spectroscopic data 

The major parameters of core-collapse supernova (SN) are (Utrobin |2003. For the above reason, hydrodynamic simula- 

g . thought to be linked to the initial stellar mass on the mam se- tions of SNe jjp have been per f orme d for only a handful of 

quence, the progenitor mass. However, the genealogy of differ- eyents: SN SN l999tm ^ SN 200 3Z, and SN 2005cs. For 

ent varieties of SNe is as yet poorly known. Fortunately, hy- these particular caseSj the inferred proge nitor masses have been 

drodynamic modeling of the light curves and the expansion ve- found unexpectedly t0 be in me range of 15 _ 22 Mq (Utrobin 

locities allows us to estimate SN parameters such as a pre-SN & chugai weU above me median yalue of K u Mq 

radius, an ejecta mass, an explosion energy, and a radioactive for the Salpeter Mtial mass distribution in me mass range of 

56 Ni amount. In the case of SNe IIP, the mass lost prior to the 9 _ 25 Mq responsible for SNe IIP (Heger et al. |2003), i.e., the 

pre-SN stage is relatively small for stars with an initial mass progenitor masses ^ on average more massive man expected, 
less than 25 Mq, so the ejecta mass provides us with a reli- 
able estimate of the progenitor mass or at least its lower limit. A more challenging problem is one related to the analy- 

Compared to other core-collapse SNe Ib/c and SNe Iln, we are sis of the sub-luminous SN 2005cs. The progenitor mass of 

able to recover the ejecta mass of SNe IIP from hydrodynamic ~ 1 8 M inf erred from hydrodynamic modeling (Utrobin & 

modeling with greater confidence because of both an accurate Chugai 12008) was found to significantly exceed the progeni- 

estimation of the photospheric velocity related to the high opac- tor mass of 6 - 13 Mq reco vered f rom archival HST images 

ity of the hydrogen-rich matter, and a low contribution of the (Maund et al. |20051 Li et al. |2006t Eldridge et al. |2007|l. No 

circumstellar interaction to the SN luminosity. reasonable explanation of this disparity has been proposed. 

Hydrodynamic modeling should only be applied to well- To pinpoint the cause of this disagreement between the 

observed SNe IIP. An adequate simulation of the bolometric two methods of the mass determination, one needs to confirm 

this discrepancy in mass measurements for a larger sample of 

Send offprint requests to: V. Utrobin, e-mail: utrobin@itep.ru SNe IIP with a broad range of observational characteristics. 
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In this respect, the well-observed luminous SN 2004et in the 
nearby galaxy NGC 6946 is a particularly favorable case. This 
object discovered soon after its explosion has the highest in- 
trinsic luminosity among well-studied events (Sahu et al. 2006 ) 
and, perhaps, the highest ejecta mass and explosion energy. In 
this case, the progenitor was directly identified in the archival 
images (Li et al. 120051 1. 

Here we perform hydrodynamic modeling of SN 2004et to 
recover the basic parameters: pre-SN radius, ejecta mass, ex- 
plosion energy, and radioactive 56 Ni mass. A brief description 
of the hydrodynamic model is given in Sect. 12.11 and the basic 
parameters of the optimal model are obtained in Sect. 12.21 In 
Sect. 12.31 we investigate whether the non-evolutionary model 
should be used instead of evolutionary pre-SN for the one- 
dimensional hydrodynamic modeling of SN 2004et and in gen- 
eral SNe IIP. The progenitor mass of SN 2004et is evaluated 
and compared to estimations for other SNe IIP (Sect. l2~4l >. The 
measured progenitor mass noticeably exceeds the mass esti- 
mated from the pre-explosion images, and this disagreement is 
discussed in Sect. 13. II In particular, we propose that the explo- 
sion asymmetry could be responsible for the disagreement in 
mass estimates and explore signatures of the explosion asym- 
metry in the SN 2004et nebular spectra (Sect. [3~2l . Finally, in 
Sect.Hl we summarize the results obtained. 

We adopt a distance to NGC 6946 of 5.5 Mpc and a red- 
dening E(B - V) = 0.41 as measured by Li et al. (120051 1. an 
explosion date on September 22.0 UT (JD 2453270.5), and a 
recession velocity to the host galaxy of 45 kms -1 following 
Sahu et al. (120061 



2. Hydrodynamic model and progenitor mass 

2.1. Model overview 

The modeling of the SN explosion is performed using the 
spherically-symmetric hydrodynamic code with one-group ra- 
diation transfer (Utrobin 2004] I2007I I, which has been applied 
previously to other SNe IIP. Utrobin d2007l > found that both this 
one-group approach and the multi-group approach of Baklanov 
et al. (12005b measured similar ejecta mass and explosion en- 
ergy for SN 1999em. The basic equations and details of the 
input physics, including calculations of mean opacities, are 
described in Utrobin (2004 ). The present version of the code 
includes additionally Compton cooling and heating. The ex- 
plosion energy is modeled by placing the supersonic piston 
close to the outer edge of the 1.6 Mq central core, which is 
removed from the computational mass domain and assumed 
to collapse to become a neutron star. The principal limitation 
of the code is that the explosion asymmetry and the Rayleigh- 
Taylor mixing between the helium core and hydrogen envelope 
(Miiller et al. |1991| l cannot be correctly treated by the one- 
dimensional model. We, therefore, study a "non-evolutionary" 
pre-SN, which takes into account the result of the mixing dur- 
ing the explosion and the shock propagation in the evolutionary 
pre-SN. The distinctive feature of the non-evolutionary model 
is a smoothed density and composition jumps between the he- 
lium core and the hydrogen envelope. 
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Fig. 1. Density distribution as a function of interior mass a) and 
radius b) for the optimal pre-SN model. The central core of 1.6 
Mq is omitted. 



The resultant structure of the non-evolutionary pre-SN in 
our optimal model is shown in Fig. Q] The pre-SN model 
is defined to be a red supergiant (RSG) with a radius of 
1500 Rq, three times larger than in the case of the normal type 
IIP SN 1999em (Utrobin [2007>. This large pre-SN radius for 
SN 2004et is implied by the broad initial peak of the bolomet- 
ric light curve shown by Sahu et al. (120061 1. The adopted mixing 
between the helium core and hydrogen envelope in the optimal 
model is shown in Fig. [2] The degree of mixing determines the 
light curve at the end of the plateau (Utrobin et al. 2007 ). The 
unmixed helium-core mass adopted for SN 2004et is 8.1 Mq, 
which corresponds to the final helium core of a main-sequence 
star of « 25 Mq (Hirschi et al. 120041 ). We note that the model 
light curve is not sensitive to any variation in the helium-core 
mass of the mixed model (Utrobin et al. 2007 ). 

2.2. Basic parameters 

The search for the best-fit model parameters is performed by 
computations of an extensive grid of hydrodynamic models. 
The optimal model should reproduce simultaneously the bolo- 
metric light curve and the photospheric velocity evolution in 
the best way. As a result, we derive the ejecta mass M em = 
22.9 + 1 M Q , the explosion energy E = (2.3 + 0.3) x 10 51 erg, 
the pre-SN radius R = 1500 ± 140 Rq, and the 56 Ni mass 
M Ni = 0.068 + 0.009 Mq. The uncertainties in the basic pa- 
rameters are calculated by assuming relative errors in the input 
observational data: 1 1% in the distance, 7% in the dust absorp- 
tion, 5% in the photospheric velocity, and 2% in the plateau 
duration. In general, the errors of derived parameters should be 
somewhat larger because of model systematic errors. However, 
the latter cannot be confidently estimated unless a more ad- 
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vanced and correct model is developed. The model uncertain- 
ties will be discussed below in Sect. 13.21 

The density distribution in the freely expanding SN enve- 
lope is shown in Fig. [3] Multiple shells in the outermost layers 
with velocities v > 14000 kmr 1 (Fig. [3} form at the shock 
breakout stage by the radiative acceleration in the optically thin 
regime. The origin of these shells is related to the specific be- 
havior of the line opacity in the outer rarefied layers of tempera- 
ture ~ 10 5 K. The innermost shell of mass ~ 5 x 10~ 3 Mq at the 
velocity of 1 1 700 km s _1 is the thin shell formed by the effect 
of the shock breakout in the optically thick regime (Grassberg 
et al. [19711 Chevalier [T98TT) . 

The optimal model describes the bolometric light curve 
quite well, including its initial (t < 30 days) peak (Fig. [4j. 
This peak is substantially broader and more luminous than the 
initial peak of SN 1999em (cf. Sahu et al. 120061 1. It is the initial 
luminosity peak of SN 2004et that requires the larger radius of 
the pre-SN model compared to the pre-SN radius of 500 Rq in 
the case of SN 1999em (Utrobin|ljo7]i. 

With one-group radiation transfer, the hydrodynamic model 
is not assumed to reproduce the monochromatic light curves 
in detail. However, it is instructive to compare the model and 
the observations in B, V, and R bands (Figs. |5^,b,c). The B 
light curve is reproduced at the initial hot phase, but not at the 
late phase. The overall fit of the V light curve is much tighter. 
In R band, the calculated light curve reproduces the plateau 
data, but does not describe the very initial stages of the light 
curve. The differences between the hydrodynamic model and 
the observations are related to deviations of the SN spectrum 
from a blackbody. These deviations are significant in the blue 
and ultraviolet at late photospheric epochs, which explains why 
the disagreement is strongest in B band. 

The computed photospheric velocity is shown together with 
two sets of observational data (Fig. [5ji): the first is recovered 
from an absorption minimum of the Fell 5169 A line (Sahu 
et al. 120061 ) and the second, from our modeling of the Na I 
doublet profile. The latter photospheric velocities can be mea- 
sured more confidently than those recovered from absorption 
minima, especially at late photospheric epochs when absorp- 
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Fig. 2. The mass fraction of hydrogen (solid line), helium (long 
dashed line), CNO elements (short dashed line), and Fe-peak 
elements including radioactive 56 Ni (dotted line) in the ejecta 
of the optimal model. 
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Fig. 3. The density and the 56 Ni mass fraction as a function of 
velocity for the optimal model at t = 50 days. Dash-dotted line 
is the density distribution fit p oc v~ 7 - 6 . 

tion lines become strong. The model photospheric velocity is 
consistent with both the Nal data and the early data of the 
Fe II 5 169 A absorption minimum. We also modeled profiles of 
the Fe II 4924, 5018, 5169 A lines at late photospheric stages, 
and found that the photospheric velocities obtained from these 
lines were rather similar to velocities obtained by modeling 
the Na I doublet profile. Unfortunately, the spectral data of SN 
2004et for the early stages are missing. The blue edge of the Ha 
absorption in the earliest spectrum on day 24 implies a maximal 
expansion velocity of ~ 12 000 - 13 000 km s _1 in the ejecta. 
This velocity is consistent with the model maximal velocity of 
12000 km s" 1 . 

2.3. Explosion of evolutionary presupernova 

The arguments in favor of the non-evolutionary pre-SN model 
leave some doubts and raise a question: why should we not 




t-ob, (days) 

Fig. 4. The calculated bolometric light curve of the optimal 
model (solid line) overplotted on the bolometric data of SN 
2004et evaluated from the UBVRI observations of Sahu et al. 
(2006) (open circles). 
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Fig. 5. Optimal hydrodynamic model. Panels a), b), and c): 
the calculated B, V, and R light curves (solid lines) compared 
with the corresponding observations of SN 2004et obtained by 
Sahu et al. (2006 ) [open circles). Panel d): the calculated pho- 
tospheric velocity (solid line) is compared with photospheric 
velocities estimated from the absorption minimum of the Fe II 
5169 A line (open circles) by Sahu et al. ((2006) and recovered 
from the Nal doublet profile (filled circles). 

consider an evolutionary pre-SN? This issue has already been 
explored for the sub-luminous type IIP SN 2005cs (Utrobin & 
Chugai 2008), for which we found that the evolutionary pre-SN 
did not allow us to produce a realistic description of both the 
light curve and the photospheric velocities. A similar problem 
was encountered while modeling the explosion of evolutionary 
pre-SNe with other hydrodynamic codes (Chieffi et al. 2003 ; 
Woosley & Heger|20Ell. 

We therefore check whether the same problem holds for 
SN 2004et, which differs in terms of both ejecta mass and ex- 
plosion energy from SN 2005cs. We adopt the pre-SN model 
with an envelope mass of 15.9 Mq and a density distribution 
that closely resembles that of the evolutionary model (Utrobin 
& Chugai 2008). The hydrogen and helium are assumed to be 




(days) t ob , (days) 

Fig. 6. Hydrodynamic model of SN 2004et for the evolution- 
ary pre-SN. Panel a): the model bolometric light curve (solid 
line) overplotted on the observational data (see Fig. |4] legend 
for details). Panel b): the model photospheric velocity (solid 
line) and the observational photospheric velocities (see Fig.[5Jl 
legend for details). 



mixed along the mass coordinate in a similar way to the op- 
timal model (Fig. An optimal fit to the bolometric light 
curve and the evolution in photospheric velocity is attained for 
an explosion energy of 1.3 x 10 51 erg and a pre-SN radius of 
600 Rq. Apparent disadvantages of the obtained model are an 
extremely narrow initial peak of luminosity and a two step- 
like transition from the plateau to the radioactive tail (Fig. |6^). 
The latter behavior is similar to that demonstrated by the light 
curves computed for the evolutionary pre-SN by Woosley & 
Heger (120071 1. The narrow initial peak is related to the rela- 
tively small pre-SN radius. However, one cannot increase the 
pre-SN radius to obtain a superior fit because the photospheric 
velocity would then become unacceptably low. A larger ini- 
tial radius would produce a higher luminosity, which, in turn, 
should be compensated by a decrease in the explosion energy, 
consequently, leading to lower expansion velocities. Even in 
the demonstrated model, the maximal velocity is only 9500 
kms 1 (Fig. |6j?), significantly lower than the observed max- 
imal velocity of 12 000 - 13 000 kms~'. We also computed 
the same model but without mixing between the helium core 
and the hydrogen-rich envelope. This model provides an even 
poorer fit because the "bump" at the end of the plateau becomes 
more boxy, in sharp contrast to the observational light curve. 

To summarize, the model including an explosion of the evo- 
lutionary pre-SN does not allow us to achieve a close fit to the 
bolometric light curve and the maximal expansion velocities 
of SN 2004et. This problem is not related to the evolutionary 
pre-SN itself. The one-dimensional hydrodynamics cannot re- 
produce the outcome of a real explosion in the evolutionary 
model, because multi-dimensional effects, in particular mixing 
between the helium core and the hydrogen envelope, play a cru- 
cial role during the explosion and shock propagation phases. 
The two-dimensional hydrodynamic model predicts that the 
Rayleigh-Taylor mixing at the helium/hydrogen interface re- 
duces the high composition and density gradients (Miiller et 
al.[T991). Our non-evolutionary pre-SN qualitatively takes this 
multi-dimensional effect into account, which ensures that the 
non-evolutionary model describes the light curve shape more 
successfully. A non-evolutionary pre-SN is also preferred by 
Baklanov et al. (120051 1 in their modeling of SN 1999em. 

2.4. Progenitor mass 

The ejecta mass M env = 22.9 + 1 Mq combined with the col- 
lapsing core of 1 .6 Mq yields the pre-SN mass of 24.5 + 1 Mq. 
A progenitor mass on the main sequence is higher by the 
amount of matter lost via the wind at the main-sequence and 
RSG phases. For the main-sequence stage, we rely on the com- 
putations of Meynet & Maeder (2003) for non-rotating stars 
with the mass-loss rate of Vink et al. (2001). They found that 
a star with an initial mass Mzams = 25 Mq lost 0.8 Mq 
during the main sequence (Meynet & Maeder 2003). For the 
RSG stage, Meynet & Maeder (2003) used the mass-loss rate 
of de Jager (119881 1 and predicted that 25 Mq and 40 Mq main- 
sequence stars lost 7.5 Mq and 21 Mq, respectively. With these 
estimates of the lost mass, we come to the progenitor mass of 
SN 2004et in the range of 30 - 40 Mq. 
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Table 1. Hydrodynamic models of type IIP supernovae. 
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A less massive progenitor is expected if we use the wind 
density recovered for the SN 2004et pre-SN from X-ray ob- 
servations. These data suggest the mass-loss rate of (2 - 2.5) x 
10~ 6 Mq yr _1 , assuming the wind velocity of 10 km s _1 (Rho et 
al. 120071 Misra et al. l2007l . Using the RSG life-time of 7 x 10 5 
yr for the 25 Mq main-sequence star (Hirshi et al. 2004), we 
find the mass lost at the RSG phase to be ~ 1 .6 Mq with an 
uncertainty of ~ +1 Mq. With the mass of 0.8 Mq lost at the 
main-sequence phase (Hirshi et al. 120041 ). the total lost mass is 
then 2.4 + 1 Mq. The pre-SN mass of 24.5 Mq combined with 
the lost mass results in the progenitor mass of 27 ±2 Mq, where 
the error includes the uncertainties in the ejecta mass and the 
mass-loss rate. The progenitor mass of SN 2004et turns out to 
be close to the maximal initial mass for SNe IIP according to 
the present-day paradigm (Heger et al. 120031) . 

Table [T] presents the parameters of all the SNe IIP studied 
hydrodynamically. The listed parameters are the pre-SN radius, 
the ejecta mass, the explosion energy, the total 56 Ni mass, the 
maximal velocity of 56 Ni mixing zone, and the minimal ve- 
locity of the hydrogen-rich envelope. The ejecta and progeni- 
tor masses of SN 2004et are found to be maximal among the 
well-studied SNe IIP. With the exception of the initial radius 
for SN 1987 A, the explosion energy and the total 56 Ni mass 
show the most extreme variations (of one order of magnitude). 
All SNe IIP are characterized by a deep mixing of hydrogen, 
indicated by the low value of v™", which is consistent with two- 
dimensional simulations (Miiller et al. I199U Kifonidis et al. 
120031 120061) . The position of SN 2004et on the plots of ex- 
plosion energy versus progenitor mass (Fig. |7h) and the total 
56 Ni mass versus the progenitor mass (Fig. strengthens the 
correlations recovered earlier for the SNe IIP studied hydrody- 
namically (Utrobin & Chugai 2008). We note that these corre- 
lations also infer the correlation between the explosion energy 
and the total 56 Ni mass, which was found and discussed earlier 
by Nadyozhin ([20031 . 

3. Discussion 

3.1. Whether the hydrodynamic mass is 
overestimated? 

The pre-SN mass recovered from the hydrodynamic modeling 
is very much close to the progenitor mass. It is, therefore, rea- 
sonable to refer to the progenitor mass thus determined as the 
"hydrodynamic" mass. Surprisingly, the hydrodynamic mass of 
the SN 2004et progenitor, 27+2 Mq, noticeably exceeds, by a 
factor of ~ 2 — 3, the value of 9+j Mq recovered from the anal- 



51.6 



op 50.6 



50.4 



-1.0 



-1.5 



-2.0 



i — 5 — i 



99em 



03Z 



I S 1 



87A 



04et 



15 



20 

MS (Mo) 



25 



30 



Fig. 7. Explosion energy a) and 56 Ni mass b) versus hydrody- 
namic progenitor mass for five core-collapse SNe. 



ysis of pre-explosion images (Smartt et al. 2009). A similarly 
large mass disparity has been found for SN 2005cs: 18 Mq 
(Utrobin & Chugai |2008l ) versus 6 - 13 Mq (Maund et al. 120051 
Li et al. 2006; Eldridge et al. 120071 ). These two cases of huge 
discrepancies in progenitor mass clearly illustrate the high level 
of uncertainty in the mass problem. 

The other side of this problem is that all the known hy- 
drodynamic masses of SNe IIP progenitors are in the range 
Mzams > 15 Mq (Fig. [7]). This mass distribution conflicts with 
the paradigm that SNe IIP originate from the mass range of 
9-25 Mq (Heger et al. 120031 ). Assuming the Salpeter initial 
mass function for SNe IIP rate in the 9-25 Mq mass range and 
neglecting selection effects, we expect that five SNe should be 
found in the 15-25 Mq mass range with a probability of only 
0.004. 

A quite different conclusion was reached by Smartt et al. 
(2009). They conclude that SNe IIP progenitors detected or un- 
detected in pre-explosion images have masses in the range of 
8-17 Mq. This significant difference in the progenitor-mass 
distributions obtained by two methods and the mass discrep- 
ancy for SN 2004et and SN 2005cs again imply that the de- 
termination of the progenitor mass remains a difficult problem. 
Although the low masses of SNe IIP progenitors from archival 
images also pose a problem for the fate of massive RSG stars 
in the range of 17 - 25 Mq (Smartt et al. 120091 ), we address 
here only the possibility that hydrodynamic masses are strongly 
overestimated compared to the real progenitor masses. 

3.2. Could explosion asymmetry be a crucial missing 
factor? 

To address the disagreements between progenitor mass deter- 
minations, we should consider possible problems with our nu- 
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merical modeling. Among the missing factors that might af- 
fect the accuracy of the inferred SN parameters, the most ap- 
parent are: multi-group radiation transfer, full non-LTE treat- 
ment of gas excitation, time-dependent ionization, high spa- 
cial resolution of the shock front, multi-dimensional effects of 
the explosion and shock propagation, consideration of density 
perturbations related to vigorous convection in the RSG en- 
velope. The effects produced by these factors perhaps differ 
in magnitude and some of the factors could be insignificant. 
However, detailed numerical studies based on advanced hydro- 
dynamic models are needed to estimate the role of each factor. 
At present, we can only state firmly that the multi-group treat- 
ment of radiation transfer cannot notably change the inferred 
SN parameters. This conclusion is based on the comparison 
between the multi-group (Baklanov et al. 120051) and one-group 
(Utrobin 2007) approaches to the SN 1999em modeling. 

A major drawback of our model may be its one- 
dimensional approximation. The multi-dimensional effects re- 
lated to the Rayleigh-Taylor mixing between the helium core 
and the hydrogen envelope during the shock propagation smear 
the composition and density jumps (Miiller et al. 119911 ). These 
effects are included artificially into our pre-SN model. More 
careful treatment of these effects with multi-dimensional radia- 
tion hydrodynamics could modify the inferred SN parameters. 

Another multi-dimensional effect, which could potentially 
be of importance, is the explosion asphericity. A growing 
amount of observational data favor a picture in which the explo- 
sion of SNe IIP is initiated by bipolar jets. SN 1987A provided 
us with the first evidence of explosion asymmetry inferred from 
polarization (cf. Jeffery 119911 ). line asymmetry (Haas et al. 
119901 ). and direct imaging (Wang et al. 120021) . During the last 
decade polarization has been detected in another five SNe IIP 
(Leonard & Filippenko 2001; Leonard et al. 2001 1 Leonard et 
al. 120061) . Two of these SNe IIP, SN 1999em and SN 2004dj, 
exhibit pronounced asymmetry in their Ha emission at the neb- 
ular stage, which is interpreted to be caused by asymmetric jets 
of 56 Ni (Chugai 120071) . The absence of a pronounced polariza- 
tion at the early photospheric epoch indicates that the explosion 
asymmetry does not lead to the asphericity of the hydrogen en- 
velope (Leonard et al. l2001|[2006l) . The spherization would de- 
velop more successfully, if bipolar jets are thermal energy dom- 
inated (Couch et al. 2009). However, the spherization of jets in 
the hydrogen envelope does not preclude that the asymmetric 
explosion could result in the modification of the velocity-mass 
distribution compared to the spherical explosion. The disagree- 
ment between the mass measurements may then be resolved, if 
the asymmetric explosion reproduces the observed light curve 
and expansion velocities for the essentially lower ejecta mass 



Table 2. Components of oxygen line decomposition. 



Component 


V 


u 


Ai 


A 2 


n 




(kms 


-') 






CUT 3 


symmetric 





1800 


0.59 


0.27 


1.25 x 10 9 


red 


1080 


900 


0.14 


0.049 


2.0 x 10 8 


blue 


-450 


480 


0.32 


0.20 


2.8 x 10 9 



£ 0.6 

(0 



1 


1 1 1 


1 1 1 1 1 1 




/ 

/ /*' 
If 

1/ 










1 


1 , 1 


1,1,1 



_l i I i I i I i I i L_ 

-4 -2 2 4 6 

v r (10 3 km s" 1 ) 



Fig. 8. Oxygen doublet [O I] 6300, 6364 A observed on day 
301 (thin solid line). Zero radial velocity corresponds to the 
rest wavelength of 6300 A. Thick solid line is the model dou- 
blet profile with all three Gaussian components, dashed line is 
the profile without the blue component, and dotted line is the 
symmetric component. 



compared to the one-dimensional explosion model. To verify 
this possibility, we would require multi-dimensional hydrody- 
namic modeling. 

Using the spectra obtained by Sahu et al. (2006}, we 
checked whether SN 2004et exhibited signatures of the ex- 
plosion asymmetry. We found that nebular Ha and [O I] 6300, 
6364 A lines indeed exhibit asymmetry. To quantify asymme- 
try effects, each line of the oxygen doublet on day 301 was 
decomposed into three Gaussian components: symmetric, red, 
and blue. The intensities of the corresponding components in 
blue and red lines of the doublet were free parameters. Because 
the line optical depth affects the doublet ratio, as in SN 1987 A 
(Chugai 1988 ), the line ratio permits us to recover the Sobolev 
optical depth and, therefore, the oxygen number density. The 
effect is weakly dependent on the electron temperature, which 
is assumed to be 5000 K. 

The result of the doublet synthesis is shown in Fig. [8] with 
the model parameters listed in Table [2] The table columns 
give the radial velocity shift of the Gaussian component (v), 
its Doppler width (u), the amplitudes of its blue (A\) and red 
(A2) doublet components, and the number density of the line- 
emitting oxygen determined from the doublet ratio of each 
Gaussian component. We emphasize that the decomposition 
is not unique unless we constrain ourselves by the shape and 
number of components. The adopted decomposition procedure 
leads to a minimal contribution of the asymmetric compo- 
nents. Although the red and blue components are weaker than 
the symmetric one, the asymmetry is rather pronounced. Both 
asymmetric components have comparable integrated fluxes, but 
are not identical in terms of the Doppler widths and velocity 
shifts. These results infer a bipolar structure and a deviation 
from the point symmetry of the line-emitting gas in the in- 
ner layers of the SN envelope. We note our conclusion refers 
to the line-emitting gas, which is not identical to the overall 
oxygen distribution. It may well be that the asymmetry of the 
line-emitting oxygen is related primarily to the asymmetry of 
56 Ni ejecta (Chugai 2007 ). Interestingly, a combination of the 
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bipolar structure and the deviation from the point symmetry is 
a specific feature of SN 1987 A, SN 1999em, and SN 2004dj 
(Chugai [200Tb . The case of SN 2004et thus provides further 
support to the conjecture that the explosion asymmetry is an 
ubiquitous phenomenon of SNe IIP. 

Remarkably, the oxygen number density for the symmetric 
and blue components (Table |2} is comparable, to within a fac- 
tor of unity, with the oxygen density of (0.5 - 1.4) x 10 9 cm" 3 
for the optimal model in the velocity range v < 2500 km s . In 
the latter case, we assume that the oxygen density is equal to 
the total matter density. This coincidence suggests that what- 
ever the role of asymmetry, the related effects do not strongly 
modify the density distribution of our optimal model in the in- 
ner layers with velocities v < 2000 km s 1 . Of course, this does 
not preclude that the velocity-density distribution in outer lay- 
ers could be modified significantly as a result of the aspherical 
explosion. 

4. Conclusions 

Our goal was to recover the parameters of the hydrodynamic 
model of the luminous type IIP SN 2004et. We obtained the op- 
timal parameter set: the ejecta mass M em = 22.9+1 Mq, the ex- 
plosion energy £ = (2. 3±0.3)xl0 51 erg, the pre-SN radius Rq = 
1500 ± 140 Rq, and the 56 Ni mass M Ni = 0.068 ± 0.009 Mq. 
The inferred ejecta mass and explosion energy are maximal 
among all the known SNe IIP explored by means of radiation 
hydrodynamics. The parameters of SN 2004et strengthen cor- 
relations between the explosion energy and progenitor mass, 
and between the total 56 Ni mass and progenitor mass discussed 
earlier (Utrobin & Chugai l2008b . 

The progenitor mass of SN 2004et, estimated by combining 
the pre-SN mass and the mass lost via the stellar wind, turns out 
to be significantly, by a factor of 2 - 3, higher than the main- 
sequence mass recovered from the pre-explosion images. This 
and the disagreement between mass estimates found earlier for 
SN 2005cs raise serious concern about the reliability of the pro- 
genitor mass recovered from the hydrodynamic modeling. We 
speculate that among the pitfalls of our hydrodynamic code, 
the most crucial could be the one-dimensional approximation. 
The artificial mixing between the helium core and the hydro- 
gen envelope, which we use to simulate real mixing, could be 
flawed, while explosion asphericity is completely ignored. The 
evidence of the explosion asphericity of SN 2004et is inferred 
from the nebular lines, which thus supports the view that the ex- 
plosion asymmetry is an ubiquitous phenomenon for SNe IIP. 
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